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Abstract. We present a computer algebra approach to proving identities on Bernoulli poly- 
nomials and Euler polynomials by using the extended Zeilberger's algorithm given by Chen, 
Hou and Mu. The key idea is to use the contour integral definitions of the Bernoulli and Euler 
numbers to establish recurrence relations on the integrands. Such recurrence relations have 
certain parameter free properties which lead to the required identities without computing the 
integrals. Furthermore two new identities on Bernoulli numbers are derived. 

Keywords: Bernoulli number, Euler number, Bernoulli polynomial, Euler polynomial, Zeil- 
berger's algorithm 

AMS Classification: 33F10; 11B68 

1. Introduction 

Bernoulli polynomials and Euler polynomials play fundamental roles in various branches 
of mathematics including combinatorics, number theory, special functions and analysis, see 
for example [U [T4l I21j . At first glance, the Bernoulli numbers, Euler numbers, and the 
corresponding polynomials do not seem to fall in the framework of hyper geometric identities. 
The powerful algorithm of Zeilberger [25] does not look like the right mechanism to handle 
the Bernoulli and Euler numbers or polynomials. 

However, as will be seen, the Cauchy contour integral representations of the Bernoulli 
numbers and Euler numbers make it possible to transform identities on these numbers and 
polynomials into identities on hyper geometric sums. In order to avoid the computation of 
the contour integrals, it is desirable to derive recurrence relations of the hypergeometric 
summands with certain parameter free properties. At the first trial, one finds it quite dis- 
appointing that the recurrence relations given by Zeilberger's algorithm seldom have the 
desired parameter free properties. Nevertheless, this drawback can be overcome by using 
an extended version of Zeilberger's algorithm. Paule [T7| first noticed that Zeilberger's al- 
gorithm can be extended to derive mixed recurrence relations for a hypergeometric term 
F(n, mi, rn-2, . . . , m r , k), where r > 1 and mi, mi, • • • , rn r are parameters. Chen, Hou and 
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Mu [3] further extended Paule's algorithm to the case with additional parameter free prop- 
erties. In this paradigm, many identities on Bernoulli and Euler numbers and polynomials 
can be verified. Moreover, some new identities can be discovered from the recurrences gen- 
erated by the original Zeilberger's algorithm without the consideration of the parameter free 
properties. 



2. Background 

Let us recall the background on Bernoulli and Euler numbers and polynomials. Let N = 
{0, 1, 2, • • • } and Z + = {1, 2, 3 • • • }. The well-known Bernoulli numbers and Euler numbers 
are defined by the generating functions 

^ ^ z n z , ^ „ z n 2e z 

l^ B n—r = — r and > E n — = . 

^-^ n\ e z — 1 n\ e zz + 1 

n=0 n=0 

By the Cauchy integral formula, we have the contour integral definitions of the Bernoulli 
numbers and the Euler numbers 

„ n\ f z dz , 
n\ f 2e z dz 

where the contour encloses the origin, has radius less than 2ir (to avoid the poles at ±2iri), 
and is traversed in a counterclockwise direction. Actually, as will be seen, there will be no 
need to compute the contour integrals, and one can formally treat the contour integrals as 
linear operators. The integral representation plays a crucial role in connecting the Bernoulli 
numbers and Euler numbers to hypergeometric terms. 

The Bernoulli numbers are also given by the following recursion 

tf"! 1 ) 5 ^ . ™>°> ( 2 - 3 ) 



k=0 



with Bq = 1. The Bernoulli numbers are rational and it is well known that -£>2n+i = for 
n > 1. The first few values of the Bernoulli numbers are as follows 

11 11 

Bn = 1 , Bt. = , Bo = — , Ba = , Br = . 

u ,1 2' 6' 30' 42 

For the Euler numbers, i?2n+i = for n > 0. 

The Bernoulli polynomials and Euler polynomials can be defined by the generating func- 
tions 

te x ^ 2e x * 

y^B n [x)— = - and V E n (x) — = - r —_ - . 

' n\ e l — 1 ^-^ n! e l + 1 

n=0 n=0 
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Clearly B n = B n (0) and E n = 2 n E n (^). The polynomials B n (x) and E n (x) obey the following 
relations 



B n( x )=ib(t) xn ~ kBk > ^ 

k=0 ^ ' 



fc=0 



We will need the following basic properties of B n (x) and E n {x). Lehmer [13J showed that 
the Bernoulli polynomials satisfy the relations B n (l) = (—l) n B n (0) and 

B n (l-x) = (-l) n B n (x). (2.6) 

Similarly, the Euler polynomials satisfy the relation 

E n (l -x) = (-l) n E n (x). (2.7) 

It is well known that the Bernoulli and Euler polynomials have the following binomial ex- 
pansions 

B n {x + y) = j^( n \B k (x)y n - k and E n {x + y) = £ (fj E k {x)y n ~ k . (2.8) 

fc=0 ^ ' k=0 ^ ' 

These basic properties will be needed for the computation of initial values for the recur- 
rence relations derived by our algorithm. 



3. The Algorithm 

In this section, we present an approach to proving Bernoulli number identities by using the 
extended Zeilberger's algorithm, and we will use an example to describe the four steps of our 
algorithm. The original Zeilberger's algorithm is devised to find recurrence relations of the 
summation Ylk -^( n > by solving the equation 

a {n)F(n, k) + ai(n)F(n + 1, k) H h aj(n)F(n + J,k) = G(n, k + 1) - G(n, k), 

where F(n, k) is a hypergeometric term in n and k, a,(n) are polynomials in n and are fc-free, 
G{n,k)/F{n,k) is a rational function in n and k. It is known that Zeilberger's algorithm 
can be applied to summands with parameters in order to establish multiple index recurrence 
relations, for example, see O Section 4.3.1] and [T7]. Recently, Chen, Hou and Mu [3] 
have found an extension of Zeilberger's algorithm to summations of hypergeometric terms 
J2k F( n -> m ii m 2-, ■ ■ ■ j "ir; k), where r > 1 and mi,m2, • • • , rn r are parameters. In fact, there 
are cases when the extended Zeilberger's algorithm becomes more efficient than the original 
form, see [3]. We will not give a rigorous description of the extended Zeilberger's algorithm, 
since it will become apparent when it is being used. 

For example, let us consider an identity of Gessel [JJ Lemma 7.2]. 
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Theorem 3.1. We have 



E(:>„ + ^(-ir-E(r) 

k — k — 



b. 



'm+k i 



(3.1) 



where m and n are nonnegative integers. 

Proof. To justify the above identity, we aim to find recurrence relations for both sides. If 
they agree with each other, then the equality is established by considering the initial values. 
There are three steps to compute the recurrence relations for the above summations. We will 
give detailed steps for the left hand side of (|3.ip . 

Step 1. Extract the hypergeometric sum from the Cauchy integral formula. 

Denote the left hand side of (13. ip by L(n,m). By the contour integral formula for B n , 
we have 



Step 2. Construct an extended telescoping equation with a shift on the parameter m of the 
summand C(n,m,k), and solve this equation by the extended Zeilberger's algorithm. 

Set the hypergeometric term 

F(n, m, k) = b C(n, m, k) + hC(n, m+l,k) + b 2 C(n + 1, m, k) + b 3 C(n + 1, m + 1, k), (3.2) 

where frj's are fe-free rational functions of n and m, namely, k does not appear in 6j's. More- 
over, we require that the rational functions 6,'s are independent of the variable z. 

By Gosper's algorithm, it is easy to check that C(n,m,k) = Zk+i — has no hyper- 
geometric solution for z^. Moreover, since the Bernoulli numbers are not P-recursive, one 
sees that Zeilberger's algorithm does not work in this case. Instead, we will try to solve the 
equation 




Denote the summand in the above integral by 




and let 



m 




F(n, m, k) = G(n, m, k + 1) — G(n, m, k) 
where G(n, m, k) is a hypergeometric term. By Gosper's algorithm, we get 



(3.3) 



r(Jfe) 



F(n, m, k + 1) 
F(n, m, k) 



a(k) c(k + 1) 
b(k) c(k) 



(3.4) 
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where 



a(k) = {m-k + l)(n + k + 1), 
6(fc) = + 1), 

c(/s) = 62 A; 2 + (62^ — &2 m + boz — 63m — b^)k — boz(m + 1) 
- b x z(m + 1) - b 2 (n + l)(m + 1) - 6 3 (n + l)(m + 1). 



Assume that G(n, m, k) = y(k)F(n, m, k), where y(k) is an unknown rational function of 
k. Substituting y(k)F(n,m,k) for G(n,m,k) in (|3.3j) reveals that y(k) satisfies 

r(k)y(k + l)-y(k) = l. 

Substituting the factorization (13, 4p into the above equation, and setting 

y(k)c(k) 



x(k) 



b(k - 1) ' 



then Zeilberger's algorithm reduces the problem further to that of finding polynomial solu- 
tions (see [18, Theorem 5.2.1]) of the following equation 

a(k)x(k + 1) - b(k - l)x{k) = c(k). (3.5) 

Notice that the coefficients a(k) and b(k) are independent of the unknowns 6,'s, and c(k) is a 
linear combination of bj's. One can estimate the degree of the polynomial x(k), as in Gosper's 
algorithm. In this case, x(k) is of degree 0. Assume that x{k) = Oq. Then the equation (|3.5p 
becomes 

(— ao — b<2)k 2 + [mb2 — nb 2 + mb^ — b^z — (n — m)ao — a$z + 63) k 
+ {in + l)(m + l)ao + + + &3^w + 6o-2 m + b^nm + n&2 
+ bizm + 6iz + 6 2 + "i^2 + mb 3 + 63) = 0. 

By setting the coefficient of each power of k to zero, we get a system of linear equations in 
ao and 6j's. Note that in the solution of this system, ao and 6j's may contain the variable 
z. To prevent z from appearing in bj's, we should go one step further to impose that the 
coefficient of any positive power of z in 6j's is zero. This may also lead to additional equations. 
Combining all these equations, if we can find a nonzero solution, then take this solution to 
the next step. Otherwise, we may try recurrences of higher order. In this case, we get a 
nonzero solution ao = — l,&o = 1, &i = — 1,&2 = 1>&3 = 0. Note that in general the 6j's are 
polynomials in n and m. 

Step 3. Compute the recurrence for L(n,m). 

By Step 2, the solution of ao, bo, . . . , 63 leads to the following telescoping equation 

C(n, m, k) — C(n,m + l,k) + C(n + 1, m, k) = G(n, m, k + 1) — G(n, m, k), (3.6) 

where 

G(n,m,k) = y, --77 ; m — rr- (-J.7) 

v ' ' 1 (k- l)!(m - £; + l)!z n+fc v y 
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Summing the above recurrence over k from to m + 1 , we obtain 



S(n, m) — S(n, m + 1) + S(n + 1, m) = 0. 



Substituting the above recurrence relation to the contour integral definition of B n , we find 
that L(n,m) satisfies 



By the same procedure, we see that the right hand side of (|3.ip . denoted by R(n,m), 
satisfies the same recurrence relation as L(n,m), namely, 

R(n, m) — R(n, m + 1) + R(n + 1, m) = 0. 

Step 4. Verify initial values. 

By considering the parity of B m , we see that (—l) m B m = B m unless m = 1. Therefore 
L(0, 1) = i?(0, 1) = 1/2 and L(0, m) = R(0, m) = B m for ra/1. This completes the proof. I 

It is known that the Bernoulli numbers and Euler numbers are not P-recursive, see [S]. 
Roughly speaking, this fact implies that the original Zeilberger's algorithm is not applicable 
to derive a recurrence relation of any order for summations involving Bernoulli numbers. For 
this reason, the extended Zeilberger's algorithm becomes necessary, and it also suggests that 
in the study of P-recursiveness of a polynomial sequence with parameters it is likely that 
one can get a recurrence relation with polynomial coefficients even for the sequence is not 
P-recursive, as long as one allows shifts on the parameters. 

4. Bernoulli Number Identities 

In this section, we give several examples of proving identities on Bernoulli numbers by using 
the extended Zeilberger's algorithm. 

The first example is the extension of Kaneko's identity given by Momiyama [15J. It was 
proved by using a p-adic integral over Z p . The Kaneko identity is stated as follows [TO] 



where B n = (n + l)B n . 

While our approach does not directly apply to Kaneko's identity because it has no pa- 
rameters, we can deal with Momiyama's identity which reduces Kaneko's identity by setting 
m = n. 

Theorem 4.1 (Momiyama's identity). 



L(n, m) — L(n, m + 1) + L(n + 1, m) = 0. 




(4.1) 




'm+k 



(4.2) 



where m and n are integers and m + n > 0. 
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Proof. Denote the left hand side and the right hand side of (|4,2p by L(n,m) and R(n,m), 
respectively. By the contour integral definition of the Bernoulli numbers, we have 



liri J e z - 1 \ ^ y ' 
J \k=0 



Denote the summand in the above summation by F(n,m,k), that is, 

Applying the extended Zeilberger's algorithm to F(n, m, k) and assuming that the output is 
independent of z, we obtain 

F(n, m, A;) + F(n, m + 1, k) + F(n + 1, m, k) = G(n, m, k + 1) — G(n, m, A;), (4-3) 

where 



G(n, m, k) 



(-l) m (m + l)!(n + fe + l) 



(k - l)!(m + 2 - k)\z n + k 
Summing the telescoping equation (|4.3p over k from to m, we are led to the following 
recurrence relation for L(n, m) 

L(n,rn) + L(n,m + l)+L(n + l,m) = -(-l) m (n + m + 2)B n+m+1 . 
Similarly, we find that R(n, m) also satisfies 

R(n, m) + R{n, m + l)+R{n + 1, m) = (-l) n (n + m + 2)B n+m+l . 
Considering the parity of B n , it is easy to see that 

((_!)- + (_i)«) ( n + m + 2)B n+m+1 = 0. 

Therefore, both sides of Momiyama's identity (|4.2p satisfy the same recurrence relation. 



To compute the initial values, setting m = we get L(n, 0) = (n + 1)-B n - It follows from 
the recursion f|2.3|) that 

E("r) B -E(:)^EG!>-o. 

fc=0 v 7 fc=0 v 7 k=0 v 7 

On the other hand, for n / 1, we have 

R( n ,o) = -(-irjr( n+ k 1 )( k + i ) B * 

k=0 V 7 
= (-l) n (n+l)B„ = {n+l)B n . 

It is easily checked that L(l, 0) = i?(l,0) = —1. So we deduce that L(ra,0) = 0) for all 
n > 0. This completes the proof. I 

The following identity is due to Gessel and Viennot [8]. 
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Theorem 4.2 (Gessel-Viennot). 

L(fc-i)/2j 1 / 2 fc_2A/2n + l\ 2n + l /2fc - 2n + 1\ 

g F^U + l JU + J g2n ' 2j = 2fc-2n + l ( v k + 1 J' n<L (44) 

Proof. Denote the left hand side and the right hand side of the above identity by L(n, k) 
and R(n,k), respectively. So we get 

L{n , k) _ J_ / _JL_ f K yf JL ( » - « )(*> + n ^-2^ 



Let 



F(n,k,j) 



1 (2k-2j\ /2n + l\ (2n-2j)! 
V fc + 1 ) \2j + l) z 2n ~ 2 i 



k — j \ k+1 

Applying the extended Zeilberger's algorithm, we get the following recurrence 

2(n + l)(2n + 3)F(n, fc, j) + 2(/s + 2)(2Jfc + 3)F(n + l,k + l, j) 

- (k + 2)(k + 3)F(n + 1, fc + 2, j) = G(n, k, j + 1) - G(n, fc, j), 

where 

G{n,k,j)- + *)!(»- « + 



(fc - 2j + l)!(2j)!(jfc + l)!z 2 ™- 2 J'+ 2 ' 

By summing the above telescoping equation over j, we obtain the following recurrence relation 
for L(n, k) 

2(n + l)(2n+3)L(n,k) + 2(k + 2)(2k+3)L(n+l,k + l)-(k+2)(k+3)L(n+l,k + 2) = 0. (4.5) 

It is easy to check that R(n, k) also satisfies the above recurrence relation. Since n < k, we 
can define L(n, n) = R(n, n) = for It is also easy to verify the initial conditions 

i f2k 4- 1 
L(0, k) = R(0, k) ' ~ ' 



2k + l\k + l 

This completes the proof. I 

It should be noted that the recurrence relation (|4.5p for L(n, k) was derived by Jacobi [9j 
in 1834, see Gessel and Viennot [8j. 

The next identity is due to Gelfand [6]. 

Theorem 4.3. We have 

B m+k mini 



~> ^\k-ln + k { ' ^\k-l 

k=i v 7 fe=i v 7 



m + k {m + n + 1)! ' 

provided that the integers m,n > are noi 6oi/i zero. 



(4.6) 
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Proof. Denote the left and right hand sides of the above identity (|4.6p by L(n,m) = 
S(n,m) +T(n,m) and R(n,m), respectively, where S(n,m) and T(n,m) are the first and 
second sums of L(n,m). Note that 

S(n,m) = — 6 — r ( > 7-irM "\ ) J^L K 



2m / e* - 1 \^ \*- V 



(n + &)z n + fe 



T ( njm ) = —<h — { > 7-ir -1 ( , n \ ^ m t k) \, )dz. 

2mje z -l\f-^ ' \k-lj (m + k)z m + k J 

Denote the summand in S(n,m) by F(n,m,k), and by the extended Zeilberger's algorithm, 
we obtain 

F(n, m, k) — F(n, m + 1) — F(n + 1, m) = G(n, m, k + 1) — G(n, m, fe), 

where 

„, , . , lNn _! m\(n + k — l)\ 

Gin, m, k) = (-1) S rrr tV; — rr- 

v ' ' 1 v y (k - 2)\{m + 2 - k)\z n+k 

Summing the above telescoping equation over k from 1 to m, we get a recurrence for S(n, m) 
S{n, m) - S{n, m + l)-S(n + 1, m) = (-l) n 

m + n + 1 

By the same procedure, or by the symmetric property T(n,m) = S(m,n), we find that 



T(n, m) — T(n, m + 1) — T(n + 1, m) = (— 1 



m + n + 1 

With the aid of the property -£>2n+i = for n > 1, we have 

L(n, m) - L(n, m + 1) - L(n + 1, m) = ((-l) m + (-l) n ) B ' m+n+1 = 0. 

' m + n + 1 

It is easy to verify that R(n, m) also satisfies the above recurrence relation. To check the 
initial values, we have 

Hn, 0, = - ± ( k I J f . £ (» ^ ) Bfc = -i, . «(„, 0). 

fc=l V 7 fe=l V 7 

This completes the proof. I 

Agoh and Dilcher [H Theorem 2.1] obtained a convolution identity for Bernoulli numbers. 
By the extended Zeilberger's algorithm and Woodcock's identity (|4,10p . we can give a direct 
proof of this result which is restated in the following equivalent form. 

Theorem 4.4. Let m, n, k > be integers, with m and k not both zero. Then 

E( n\ k\m\ , , 
. \B k+ jB m+n -j = — , , , 7T, [n + o{m, k)(m + k + l))B m+n+k 
. \J J (m + k + 1)! 



j=0 

JH { )r iWl-, ( _ 1)fc A + 1\ + _ ™] Z7 

^ v ' m + k + l-r y ' \ r )\ k + 1 k + lJ 

m+k 



n+r— 1 



Sp, x y B m+k+l-r / i)m( m + 1 \ ( m + 1 ~ r n rfc Aff 
^ J m + £; + l-/ V r A m + 1 m + 1 > 

r=0 



n+r— 1 ) 

(4.7) 
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where 6(m, k) = when m = or k = 0, and 5(m, k) = 1 otherwise. 



Proof. Let L(n, m, k) and i?(n, m, A;) denote the left hand side and the right hand side of 
the above identity (|4.7p . respectively. Our approach leads to the recurrence relation 

S(n, m + 1, k) - S(n + 1, m, k) + S(n, m, k + 1) = 0, (4.8) 

where m/0 and k ^ 0. Considering the parity of the Bernoulli numbers, we have (— l) k Bj. = 
Bk for k ^ 1. The known convolution identity on Bernoulli numbers 

(t) B k B n _ k = -nB n ^ - (n - l)B n , n > 1 (4.9) 



k=0 v 



yields that 



L(0,m,l) =B x B m = ~\s m , 



R(0, m, 1) = - ^-B rn+l + f>l) r ^f^-(-l) H ( " ^) 
m + 1 m + 1 — r \r I V 2/ 



r=0 

m+1 



+ V(-i)- Bm+2 - r (-ir r + 1 V - — U-i 

r=0 



1 m (—l) m ^-\ /m + 1\ 

—rB m+ i — B m+ i + B m B\ H — > (— l) r LB TO+ i_ r 5 r 

m + 1 m + 1 m + 1 ^-^ \ r I 

r=0 

(_ l) m /^ir 1 /m + l\ \ 
Bm+i + -Bm-Bi H — r ) ( B m+1 _ r B r - 2(m + l)B m Bi - B m+1 



r=Q 

m 



— — B m+ i + B m B\ — (—l) m B m+ i 

This gives the proof for (|4.7p when m/0 and k ^ 0. 

Moreover, if m = or k = 0, we can simplify the identity to an equivalent form of a 
known identity discovered by Woodcock [23J 

— ^2 ( u j ( — l) k B m -kB n -i+k = — ^ f ^_ J (— ^) k B n -kB m -i+k- (4-10) 
m fc=i ^ ' n fc=i ^ ' 

This completes the proof. I 



5. Bernoulli Polynomial Identities 

In this section, we show that our approach is also valid for proving identities on Bernoulli 
polynomials. We will explain how this method works by considering an identity due to Sun 

[2D]. 
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Theorem 5.1. We have 

("!)*£ (% k ~ j B l+j (y) = (-l) l Y,(% l ~ 3 B k+J (z), (5.1) 

j=0 j=0 ^ 

provided that x + y + z = 1 . 

Proof. Denote both sides of the above equation by L(k, I) and R(k, I), respectively. We have 
Let F(k,l,h,j) denote the summand in the above integral, that is, 



Applying the extended Zeilberger's algorithm to F(k,l,h,j) with the assumption that the 
output is independent of the variables u and h, we arrive at the relation 

xF(k, I, h,j) + F(k + 1, 1, h,j) + F(k, I + 1, h,j) = G(k, I, h,j + 1) - G(k, I, h,j), (5.2) 

where 

G(k,l,h,j) = —?1—F(k,l,h,j). 
k-j + l 



Summing both sides of (15. 2\\ over h and j gives the recurrence relation 

xL(k, I) + L(Jfc + 1, 1) + L(fc, Z + 1) = 0. 

Similarly, it can be shown that R{k, I) satisfies the same recurrence relation. It remains to 
check the initial values 

L(0,0 = ^(y), 

11(0,1) = (^x^B^z) = {-i)ijr( l \x l -i(B + z y 

= (-l) l (B + x + z) l = (-l) l Bt(x + z) = (-1)^(1 -V) = Bi(y), 

as desired. I 

It is worth noting that the extended Zeilberger's algorithm is indeed efficient in deriving 
recurrence relations for multiple sums. The next identity is given by Wu, Sun and Pan |24j . 

Theorem 5.2. We have 

m 



jn / _i_ -i \ 

:-irE( m r )(n+k+i)B n+k ( X ) 

+ (-1)" E ( n t *) (™ + k + l)B m+k {-x) 

= {-l) m {n + m + l){n + m + 2)x n+m . (5.3) 
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Proof. Denote the two sums on the left hand side of (|5.3p by S(n, m) and T(n, m) respec- 
tively. Let L(n, m) = S{n, m) + T(n, m), and let R(n, m) denote the right hand side of (15. 3p . 
Write 



Denote the summand in the above expression by F(n,m,k,j). Applying the extended Zeil- 
berger's algorithm with the assumption that the output is independent of the parameters z 
and j, we obtain that 

F(n, m, k,j) + F{n + 1, m, k,j) + F(n, m + 1, k, j) = G(n, m,k + 1, j) — G(n, in, k,j), 

where 

k 

G(n,m,k,j) = -——F(n,m,k,j). 

m — k + 1 

By summing the above telescoping equation over j from to n + k and k from to m + 1, 
we deduce that 

S(n, m) + 5(n + 1, m) + S(n, m + 1) = (-l) m+1 (n + m + 2) J B n+m+ i(x). (5.4) 

From the symmetry property it follows that T(n,m)(x) = S(m,n)(—x). This leads to the 
following recurrence relation for T(n, m) 

r(n,m)+T(n,m + l)+T(n + l,m) = (-l) n+1 (n + m + 2)B n+m+1 {-x). (5.5) 

Adding (|5.4p to (|5.5p . we derive a recurrence relation satisfied by L(n,m) 

L(n, m) + L(n + 1, m) + L(n, m + 1) 

= (-1)™+^ + m + 2)£ n+m+1 (x) + (-l) n+1 (n + m + 2)B n+m+1 {-x) 

n+m+l 



■»- t-m- t-i / i I 1 \ / 

(_l)m+l (n + m + 2) ( n + ^ +i )x" +m+1 - fe J B fe (l + (-l) fc+1 

n+rra+1 , i 1 \ 

2(_i)-+i( n + m + 2 ) [ k )x n+m+1 - k B k 

fe=n \ / 



= 2(-l) m+1 (n + m + 2)(n + m + l)x n+m Si 
= (-l) m (n + m + l)(n + m + 2)x n+m . 

It is easy to see that R(n,m) satisfies the same recurrence relation as L(n,m). Based on the 
well-known identity for Bernoulli polynomials 



™ n ~ l = E uH-^) = E 

fc=l V 7 fc=0 V ' 



it is straightforward to verify that 



L(n, -1) = + (-1)™ £ ( n + ^ 
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n-1 , v 
(-!)>+!)£ ^(-^) 
fc=0 ^ ' 

(-l) n (n + l)n{-x) n - 1 = -n(n + l)^" 1 = R(n, -1). 



This completes the proof. I 

Note that the above identity (15. 3p reduces to Momiyama's identity (14, 2|) by setting x = 0. 
We also note that integrating the identity (|5.3p over x and using the Bernoulli number identity 
(I3.1|) . one can derive the following identity of Wu, Sun and Pan [21] 

m / \ n / \ 

(-ir E n = E n B ™+i(-x)- (5-6) 

i=0 v y j=0 VJ/ 

The following identity is derived by Sun [20J. 
Theorem 5.3. We have 

( } h\j) i+j + i +i ) h\i) + i ~ (k + i + i)( k i l Y (5 ' 7) 



3=0 w J j=0 

provided that x + y + z = 1 . 



Proof. Let L(k,l) and R(k,l) denote the left hand side and the right hand side of (|5.7p . 
respectively. It can be shown that 

xL(fc, Z) + L(Jfe + 1, Z) + L(k, I + 1) = 0. 

It can also be shown that R(k, Z) satisfies the same recurrence relation. To check the initial 
conditions, we have 



L(0,Z) 



B i+i{y) . f f l \ j-j I3 i±ii z l 



^Tr + TTT^( v z- J + iJ a;i? ^ +l( 



j^-B l+1 (y) -tHa* 1 (by m 
Z + l Z + l + yy> Z + l y J M — 

i?(0,Z), 



Z + l 
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as desired. 



We remark that the above identity (|5.7p reduces to (|5.ip by viewing z = 1 — x — y as a 
function of y and by taking partial derivative with respect to y. It also specializes to (j5.6j) 
when setting y — > x and z = —y — * —x. Moreover, differentiating both sides of (|5.7p with 
respect to y twice, we obtain the following identity derived by Sun [20], which can be verified 
by our approach. The proof is omitted. 

Theorem 5.4. Suppose that x + y + z = 1, then 

3=0 V J 
I 



= (-l) fc (A; + I + 2)(B fc+l+1 (x + y) - B k+l+1 (y)). (5. 



In [21 j. Theorem 1.1], Sun and Pan find a symmetric relation between products of the 
Bernoulli polynomials. 

Theorem 5.5. Let n € Z + and x + y + z = 1. If r + s + t = n, then 

r B- 1 )'(0(„!t) J ^ (x) * (,) 

fc=0 v / v / 

+»g(- i ) fc Q)( n :,)^w^w 



*t(-D'(0 

fc=0 v 7 



s 

n — k 



B n _ k (z)B k (x) = 0. (5.9) 



Proof. Denote the three sums on the left hand side of the above identity by S(n,r,s), 
T(n, r, s), R(n, r, s) respectively. Since n = r + s + t, S(n,r,s) can be expressed as 

^)7?^/s^(»- 1 K3("--r)("7*)(!D5^ , '")' fa *- 

Our approach yields the following recurrence relation 

(s + l)S(n, r + 1, s) + (r + l)S(n, r, s + 1) + (n - r - s - l)5(n, r + 1, s + 1) = 0. 

Similarly, it can be shown that T(n, r, s) and i?(n, r, s) satisfy the same recurrence relation. 
Since < r, s < n and n £ Z + , we obtain that 



,S'(//. 0..<! + ■:/'( //.(). .s)+ R( /?.<).*) - S j5 n (z) + (n-s)r )/?„(-.) - U. 



and 

S(n, r, 0) + T(n, r, 0) + R(n, r, 0) = r ~ ^ B n (i) + (n - r)(-l) n Q 5 n (a?) = 0. 
It follows that S(n, r, s) + T(n, r, s) + i?(n, r, s) is identically zero. This completes the proof. 
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6. Euler Number and Polynomial Identities 



In this section, we show how to prove identities on Euler numbers and polynomials by using 
our approach. As the first example, we consider the following identity due to Wu, Sun and 
Pan pi]. 

Theorem 6.1. We have 

(-d-E = E (") E ™« H) • i"' 

where m and n are nonnegative integers. 

Proof. Denote the left and right hand sides of (|6.ip by L(n,m) and R(n,m), respectively. 
By the contour integral definition of the Euler numbers (12. 2D and the relation (|2.5p . we have 



\ fc=0 / 



, j=0 fe=0 

Denote the summands in the above two integrands by 

/ m\ {n + k)\ 



S(m, n, k) = (— 1) 



k I 2 n ~^^ z n +^+i 
n fn\fm + j\. . m+ j_ k kl 



2 fc z fc+r 

Applying the extended Zeilberger's algorithm, we obtain 

S(n, m, k) + S(n, m + 1, k) + S{n + 1, m, k) = G(n, m, k + 1) — G(n, m, k), 

T(n, m, fc, j) + T(n, m + 1, A;, j) + T(n + 1, m, k,j) = H(n, m, k,j + 1) — H(n, m, k, j), 

where 

G(n,m,k) = -S(n,m,k), H(n,m,k,j) = — -4. : T(n,m, k, j). 

m+ 1 — k n+ 1 — j 

Therefore, L(n, m) and R(n,m) satisfy the same recurrence 

L(n, m) + L(n, m + 1) + L(n + 1, m) = 0. 
Consequently, the identity (|6.ip can be verified by computing the initial values 

m -) - (-D- E (?) § - E (T) (-D-f - - 5) - 

fe=0 v 7 fc=0 v 7 

as desired. I 

Wu, Sun and Pan [M] also derived an identity by substituting the Bernoulli polynomials 
in (|5.6p with Euler polynomials. This identity can be verified by our approach. The proof is 
omitted. 
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Theorem 6.2. We have 



m / \ n / \ 

(-i) m £(fc)w*) = (-i) B £(£)^ 

fc=n ^ ' fc=n ^ ' 



? m+fe (-x). (6.2) 

k=0 v '" / fc=0 

Note that differentiating both sides of the identity (|6.2p with respect to x leads to the 
following identity also due to Wu, Sun and Pan 



H) m E I, )(n + fc + l)^(x) 
fc=o ^ ' 

l — n V / 



k=0 

n+m+1 



l)- 2 (n + m + 2)(x n+m+i - E n+m+1 {x)). (6.3) 



The following identity is derived by Sun [20 
Theorem 6.3. 



i=0 w — i=0 w — (*+i+i)m 

provided that x + y + z = 1. 



Proof. Denote the two sums in the left hand side of the above identity by S(k, I) and T(k, I) 
respectively. Let L(k, I) = S(k, I) +T(k, I), and let R(k, I) denote the right hand side of (|6,4p . 
By computation, we find 

xS(k, I) + S{k + 1,0 + S(k, I + 1) = 0. 

Since T(k, I) = S(l, k), 

xT(k, I) + T(k + 1,1) + T(k, l + l) = 0. 

Therefore, 

xL(k, I) + L(k + 1,0 + L(k, I + 1) = 0. 

It is easy to check that R(k, I) satisfies the same recurrence relation. To check the initial 
values, we have 

^+m>^£(T> ,+i ^ (z) 

7 = 1 



Ei. 



— + (-l) l T^-(E l+1 (x + z)-x w 



l + l K ' 1 + 1 
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--^ -12(0,0, 

as desired. I 

Our approach can also be applied to identities involving products of the Euler polynomials 
and the Bernoulli polynomials. We take the following identity of Sun and Pan [211 Theorem 
1.1] as an example. 

Theorem 6.4. Let n £ Z + , r, s, t > 0, r + s + t = n — 1 and x + y + z = 1, then 

B-i>*Q( n :>*M^) 



(-«"B-«'(l)(„!^w^« 



iD-i)'(;)(„_ t 1 _ i )^(^„- 1 - i w. (6.5) 



Proof. Denote the two sums on the left hand side of the above identity (|6.5p by S(n,r,s) 
and T(n,r,s) respectively. Let L(n,r,s) = S(n,r,s) — T(n,r,s), and denote the right hand 
side of (|6.5p by R(n,r,s). Note that 



fc=0 j=0 h=0 



x x fc J ^r[ '", " I [z - - ) T—TT-, \dudv 



\ \ n—k—h [ 



v?\ h J\ 2 J 2 h v h+1 

Applying the extended Zeilberger's algorithm, we have 

(s + l)S(n, r + 1, s) + (r + l)5(n, r, s + 1) + (n - s - r - 2)S(n, r + 1, s + 1) = 0. 

It can also be shown that T(n, r, s) and i?(n, r, s) satisfy the same recurrence relation. Since 
< s < n — 1 and < r < n — 1, it follows that 

j£„(z)-(-lW ^ =0 = R(n,0,s), 

L(n,r,0) = (-1) b Qb b (i) - E("l) n+fc Q (" ~ 1 ^^(y)^^) = = fl(n,r,0). 
This completes the proof. I 



7. Deriving New Identities 

Applying the original Zeilberger's algorithm to a Bernoulli number summation, we may obtain 
a recurrence relation for the summand which contains the integral variable z. Although such 
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a recurrence cannot be used to prove the Bernoulli number identity itself, it may lead to a 
new identity. For example, let us consider Kaneko's identity (|4.ip 

n + l 



k=0 



where .B n = (n + l)i? n . From the recurrence obtained by Zeilberger's algorithm, we can get 
the following generalization of this identity. 

Theorem 7.1. VFe /icrae 

n+3 



)(n + k + 3)(n + k + 2)B n+k = 0. (7.1) 
fc=o ^ ' 



Proof. Denote the left hand side of Kaneko's identity by L(n). By the contour integral 
definition of the Bernoulli numbers, we have 



n+l / i i\ 



Denote the summation in the above integral by S(n). Obviously, 

L(n) = — I —^—S(n)dz = 
2iri J e z — 1 

for all n > 0. Applying Zeilberger's algorithm, we get 

z 2 S(n + 2) = 2(n + 3)(2n + 5)5(n + 1) + (n + 2)(n + 3)5(n). 
By integrating over z on both sides of the above recurrence, it follows that 

— I —^—z 2 S(n + 2)dz 
2-Ki J e z - 1 v ' 

f t / i (ra + A; + 3)- dz 



vri J e z - 1 \ ^ 

n+3 / m 

E( n I J (n + A; + 3)(n + k + 2)(n + /c + l).B n+ & 
fc=o ^ ' 

11 [l J (n + k + 3)(n + k + 2)B n+k 

2(n + 3)(2n + 5)— / — — STn + l)cfc + (n + 2)(n + 3)— / — ^S"(n)dz = 0. 
27n J e z — 1 2m J e z — 1 



This completes the proof. 
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Gessel Theorem 7.3] extended Kaneko's identity (|4.ip to the following form 

n+l / , n \ m—1 

n + 



— E m " +1 A T )^n+fc= E ((2n + l)A:-(n + l)m)A ; n (A ; -mr- 1 . (7.2) 
fc=o ^ ' fc=l ^ ' 



Note that when m = 1, the above identity becomes Kaneko's identity. From the above 
identity, we can deduce the following theorem by applying Zeilberger's algorithm. 

Theorem 7.2. We have 

n+3 / , \ m—1 

. ,,, / . ••• \ i i • -,^,,-t-i. / ,/ '•• \ *" m) n , (7-3) 

(n + 3) ' 



J- m ^- k r + j {n + k + 3) (n + k + 2)B n+k = ^ p(n, m, k)k n (k 
' k=o ^ ' k=l 



where 

p(n, m, k) =2(n + 2)(2n + 3)(2n + 5)k 3 - 2m(n + 2)(2n + 5)(3n + b)k 2 
+ 3m 2 (n + 2)(2n 2 + 7n + 7)/c - m 3 (n + l) 2 (n + 2). 

Proof. Denote the left hand side and the right hand side of (I7,2p by L(n,m) and R(n,m), 
respectively. Then we have 



L(n, m) = J2 ™ n+1 ~ k ( n+ k l \ n + k + l ) B n+l 



Im J e z - 1 \ ^ 

J \ k=0 



■ m n+i-k (n + A (n + fc + l) (n + k)\ \ 
2m J e z - 1 \ ^ I k J n + l z n+k ) 



Denote the summation in the above integral by S(n,m). By Zeilberger's algorithm, we find 
that 

z 2 S{n + 2, m) = 2(n + 2)(2n + 5)5(n + 1, m) + m 2 (n + l)(n + 2)5(n, to). (7.4) 
Integrating the left hand side of the above recurrence over z, we get 

— / — - — (z 2 S(n + 2, m))dz 
2m J e z - 1 v v ' ;; 



+3-fc p + 3\ (n + ft + 3) (n + 2 + fc)! \ & 
/ n + 3 z n+fc I 



2vri / e z - 1 I ^ 

J \ k=0 

1 n+3 / _i_ o\ 

= r^E mn+3 ^l I )(n + k + 3){n + k + 2)B n+k . 

On the other hand, integrating the right hand side of (|7.4p over z and substituting L(n, to) 
by R(n,m), we obtain 



— / — - — ( 2(n + 2)(2n + 5)5(n + 1, to) + m 2 (n + l)(n + 2)S(n, m) ] cte 
2-7ri / — 1 \ / 

2(n + 2)(2n + 5)L(n + 1, to) + m 2 (n + l)(n + 2)L(n, m) 
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2(n + 2)(2n + 5) ^ I (2ra + 3)fc - (n + 2)m ) k n+l (k - m) 




k=l 



n 



+ m 2 (n + l)(n + 2) ^ I (2n + l)k - (n + l)m ) A; n (A; - m) 




fc=i 



n-1 



m—1 



^ m, k)k n (k — m) 

k=l 



n-1 



as desired. 



Obviously, the above identity (17,3j) reduces to (|7.1j) by setting m = 1. 

8. Concluding Remarks 

To conclude this paper, we remark that our approach is not restricted to identities on Bernoulli 
and Euler polynomials. It also applies to sequences ao,a%,a2, ■ ■ ■ whose generating functions 
f(z) lead to contour integral representations of a n with hypergeometric integrands. For exam- 
ple, the Genocchi numbers fall into this framework. We can apply the extended Zeilberger's 
algorithm to prove the following identity on Genocchi numbers 



where m, n £ Z + . Recall that the Genocchi numbers can be defined by the generating 
function 



We note that there are other approaches to proving identities related to special numbers 
and functions. For example, Paule and Schneider |16j used Karr's summation algorithm in 
difference fields and Zeilberger's algorithm to prove harmonic number identities and de- 
rive new identities. Kauers [12j gave an algorithm which can be applied to verify many known 
identities on Stirling numbers and to discover new identities. Stan [19] applied Wegschaider's 
mathematica software package MultiSum [22] to deal with identities related to Poisson in- 
tegrals. Moreover, the package MultiSum can establish multiple index recurrence relations 
for the hypergeometric terms with parameters which can also be established by using the 
extended Zeilberger's algorithm. 
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n / \ m s \ n+k 

E(;)(-i)'GU*-E(T)(-i)*D-i) , G. 

k=0 V 7 k=0 V 7 i=0 



n=l 
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